STABILIZATION OF THE WAVE EQUATION WITH EXTERNAL 

FORCE. 



M. DAOULATLI 

Abstract. We study the rate of decay of the energy functional of solutions of the wave 
equation with localized damping and a external force. We prove that the decay rates of the 
£Nj ' energy functional is determined from a forced differential equation. 
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This article is devoted to the study of stabilization for the wave equation with external 
force on a compact Riemannian manifold with boundary. In the first part of this paper, we 
consider the following wave equation with linear internal damping and external force 

d?u - An + a (x) d t u = f (t, x) R + x M 

u = R + xdM (1.1) 

(u(0),3tu(0)) = (uo,«i) 



£f) ' Here M = (M, q) is a compact, connected Riemannian manifold of dimension d, with C°° 

boundary dM, where q denotes a Riemannian metric of class C°°. A the Laplace-Beltrami op- 
erator on M. a (x) is a non negative function in C°° (M) and / is a function in L 2 (R+ x M) . 

00 ' We define the energy space 

n = Hi (M) x L 2 (M) 
where 

Hi (M) = {ueH 1 (M) ■ u\ dM = 0} 

which is a Hilbert space. Linear semigroup theory applied to (jl.ip . provides the existence of 
a unique solution u in the class 

u G C° (R+, i^o 1 (M)) n C 1 (R+, L 2 (M)) 

With (jl.ip we associate the energy functional given by 

K(t) = - / |v-u(t,x)| 2 + |a t u(t,x)| 2 dx. 

^ Jm 

The energy ^(t) is topologically equivalent to the norm on the space %. Under these as- 
sumptions, the energy functional satisfies the following identity 

E u {t)+ [ fa (x) \d t u\ 2 dxda = E u (s)+ [ [ fd t udxda (1.2) 
is Jm Js Jm 

for every t > s > 0. 
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The topic of interest is rate of decay of the energy functional. This problem has a very long 
history. The connection between controllability, observability and stabilization was discovered 
|24j and effectively used in the context of linear PDE systems. 

When / = and the damping term acts on the hole manifold the problem has been 
studied by many authors [;2j and reference therein. For the wave equation with localized 
linear damping term, we mention the works of Rauch- Taylor [22] and Bardos et al [3] in 
which microlocal techniques is used. In particular the notion of geometric control. We cite 
also the works of Lasiecka et al [16] and Triggiani- Yao [26] in which another approach based 
on Remannian geometry is presented. 

Particular attention has been paid to the case when M is a bounded domain and the 
damping is linearly bounded [13] P3] and reference therein. Under certain geometric condition, 
the energy functional decays exponentially. Damping that does not satisfy such linear bound 
near the origin (e.g. when the damping has polynomial, exponential or logarithmic behavior 
near the origin) results in a weaker form of the energy that could be expressed by algebraic, 
logarithmic (or possibly slower) rates [15] [1] [19] . Finally we mention the work of cavalcanti 
et al [6] when M is a compact manifold with or without boundary. 

When / 7^ 0, the literature is less furnished, we specially mention the works of Haraux [12] 
and Zhu [27] when the damping is globally distributed. 

We should also remark when the support of the dissipation may be arbitrarily small require 
more regular initial data and result in very slow (logarithmic or slower ) decay rates as shown 
in [7] and reference therein. 

We assume that the geodesies of M have no contact of infinite order with dM. Let oj be 
an open subset of M and consider the following assumption: 

(G) (u, T) geometrically controls M, i.e. every generalized geodesic of M, travelling with 
speed 1 and issued at t = 0, enters the set w in a time t < T. 

This condition is called Geometric Control Condition (see e.g. [4]) We shall relate the 
open subset uj with the damper a by 



Under the assumption (G) it was proved in [31 EH], that the energy decays exponentially, 
moreover if there exits a maximal generalized geodesic of M that never meets the support of 
the damper a, then we don't have the exponential decay of the energy for initial data in the 
energy space. 

It is known that the exponential decay of the energy is equivalent to the following observ- 
ability inequality: 

(A) Linear Observability inequality: There exist positive constants T and a = 
a(T), such that for every initial condition ip = (uq, u±) G H the corresponding solution 



for every t > 0. 

In this paper, under the assumption (G), we show that for the non autonomous case the 
corresponding observability inequality reads as follows: 

(B) Non autonomous linear Observability inequality: There exist positive con- 
stants T and a = ot(T), such that for every initial condition p = (uq,U\) G T~L the 



{x G M : a(x) > 0}. 



satisfies 




■t+T 



(1.3) 
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corresponding solution satisfies 

ft+T r 

E v (t)<a / a{x) \d t v\ 2 + \f {s, x)\ 2 dxds. (1.4) 

Jt Jm 

for every t > 0. 

From the observability inequality above, we infer that the rate of decay of the energy will 
depends on J M \f (t, x) | 2 dx. Now we state the main result of the first part of the paper: 

Theorem 1. Let u(t) is the solution to the linear problem (COP with initial condition 
(u ,ui) G %. We assume that (w,T) satisfies the assumption (G) and 

r(t) = d, T / |/(t,x)| 2 dxGL 1 (M + ) 
Jm 

with Ci } t > 1. Then 

E u (t) < 4e T ^5 (t - T) + y T (s) dsV t > T 
where S (t) is the solution of the following ordinary differential equation 

^ + J^ = r(t), 5(0) = E„(0). (1.5) 

where Ct > 1. 

1.1. Applications for the linear case. Setting 

T(t) = C\ T [ \f(t,x)\ 2 dx 
Jm 

with Ci t T > 1. 

The ODE (|1.5p governing the energy bound reduces to 

^+CS = r(t) (1.6) 

where constant C > does not depend on E u (0) . 

(1) If there are constants M > and 6* > 0, such that 

r (t) < Me" 9i 



We have 

ft i 



e 9T -l 



h-T l 

Multiply (|1.6p both sides by exp(Ci) and integrate from to t, we obtain 

(a) C > e 

E u (t) <c(l + E u (0))e- et ,t >0 

(b) c = e 

E u it) <c(l + E u (0)) (1 + t) e' et , t > 

(c) C < 9 

E u (t) <c{l + E u (0))e- ct ,t>0 
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(2) If there are constants M > and 8 > 1, such that 

r (i) < M (1 + 

We have 

/ (i + s y e ds < T(i + t-Ty e ,t > t 

Jt-T 

In order to obtain the rate of decay in this case, we use proposition [TJ Then 

E u (t) <c(l + t-Ty e ,t>T 

where c > and depends on E u (0) . 

Remark 1. If we consider the following system 

dfu -Au + a(x)g (d t u) = f {t, x) R + x M 

u = R + xdM 
{u (0) , 8 t u (0)) = («o, u x ) € Hi (M) x L 2 (M) 

with g continuous, monotone increasing function, vanishing at the origin and linearly bounded. 
Then the result of the theorem above remains true. 

1.2. The nonlinear case. In the second part of the paper we study the rate of decay of 
the energy functional of solution of the wave equation with nonlinear damping and external 
force. More precisely, we consider the following system 

dl u - Au + a (x) g (d t u) = f (t, x) R+ x M 

u = R + xdM (1.7) 

(u (0) , dtu (0)) = (uq, u\) € tfo 1 (M) x I? (M) 

g is a continuous, monotone increasing function vanishing at the origin. Moreover we 
assume that, there exists a positive constant m, such that 

— \s\ 2 < g (s) s < m\s\ 2 , \s\ > rj (1-8) 
m 

for some r] > 0. a (x) is a non negative function in C°° (M) and / is in I? (R + x M) . 

Nonlinear semigroup theory applied to (jl.ip . provides the existence of a unique solution u 
in the class 

u € C° (R+, ^ (M)) n C 1 (R +1 L 2 (Af)) 
Under these assumptions on the behavior on the damping, the energy functional satisfies the 
following identity 

E u (t)+ / / a (x) g (d t u) d t udxdcr = E u (s) + / / fd t udxda (1.9) 

Js J M J s J M 

for every t > s > 0. 

It is well known, for the nonlinear problem without a external force the corresponding 
observability inequality |151 [9].., reads as follows: 

(C) Nonlinear Observability Inequality: There exists a constant T > and a con- 
cave, continuous, monotone increasing function h : R+ — > R_|_, /i(0) = (possibly 
dependent on T ) such that the solution u(t, x) to the nonlinear problem (jl.7p with 
initial data = (uo,u\) and / = satisfies 

£ u (t) ^ h (yj t j a{x)g(d t u)d t u dxds^j , (1.10) 
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for every t > 0. 

The function h(s) in (jl.lOp depends on the nonlinear map g(s), and ultimately determines 
the decay rates for the energy E u (t). The energy decay for the nonlinear problem will be 
determined from the following ODE 

s t + hr l (cs) = o, s(o) = e u (o) (i.ii) 

we show that under the assumption (G) we obtain the following observability inequality 

(D) Nonlinear Non-autonomous Observability Inequality: There exists a con- 
stant T > and a concave, continuous, monotone increasing function h : M+ — > R+, 
/i(0) = (possibly dependent on T ) such that the solution u(t, x) to the nonlinear 
problem (|1.7j) with initial data <p = (uo,u\) satisfies 

E u (t)<h( f [ a{x)g(d t u)d t u dxds + [ [ \f {s,x)\ 2 dxds) , (1.12) 
\Jt Jn Jt Jm J 

for every t > 0. 

Before giving the main result of this section, we will define some needed functions. Ac- 
cording to [15] there exists a strictly increasing function ho with ho (0) = such that 



h (g (s) s) > e (js| 2 + \g (s)\ 2 ^j , \s\<rj 



for some €o,tj > 0. For the construction of such function we refer the interested reader to 
[Hill]. With this function, we define 

h = I + m a (M T ) h o — — 

m a (M T ) 

where 

m a = a (x) dxdt and M T = (0, T) x M 
We can now proceed to state the main result of the second part of the paper 

Theorem 2. Let u(t) is the solution to the nonlinear problem j 1. 7| ) with initial condition 
{uo,u\) S Tl. We assume that (uj,T) satisfies the assumption (G) and 

Tit) = 2 J m \f (t, x)\ 2 dx + r(\\f (*, •) \\l*(m)) G L 1 (M+) 

where tp* is the convex conjugate of the function ip, defined by 

i _ — -i <' ° 2 



with Ct > 1. TTien 



+oo s G 



(t) < 4e T (s(t-T) + J r (s) dsV t > T 
where S (t) is the solution of the following ordinary differential equation 



d ~§ + jfh- 1 (±s) = r (t) , S (0) = E u (0) . (1.13) 



with, K >C T - 
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1.2.1. Applications for the nonlinear case. 

Proposition 1. Let p a differentiable, positive strictly increasing function on R + , We assume 
that there exists mi > such that, p (x) < mix for every x € [0, n] for some < rj « 1 and 
the following property 

p (Kx) > mp (K) p (y) (1-14) 

holds for some m > and for every (K,x) € [1, +oo[ x R + . T G C 1 (M+) • Let S satisfying 
the following differential inequality 

^+p(S)<T(t), 5(0)>0. 

(1) r (t) = /or every t > 0. We assume that 5 (0) > 0. Then 

S (t) < V" 1 (*) , /or every i > 

where 

x e 10,5(0)] . 

(2) r (t) > /or every t > 0. 

(a) There exist c > and re > 1 such that 

d -p~ l (r (t)) + cr (t) < 0, for every t > (1.15) 



dt 



and 

Then 
where 



mp(re) - rec - 1 > . . 

rep- 1 o r (o) > 5 (o) li,iDj 

5 (i) < reV'^ 1 (ct) , /or every t > 
rp-ior(o) ds 



Y> (x) = / 



p(s) 

x € jO,^ 1 oT(0)] . 
(b) There exist c > and re > 1 swc/i i/iai 

-^-p -1 (r (t)) + cr (t) > 0, for every t > 

and 

mp (re) — ere — 1 > 
re^ 1 oT (0) > 5(0) 

JTten 

5 (t) < rep" 1 o T (i) , /or every t > 

Proo/. 

(1) If 5(0) = 0, since 5 is positive and decreasing then 5 (t) = for every £ > 0. We 
assume that 5 (0) > 0. Let ip, the function defined by 
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then if) is a strictly decreasing function on (0, S (0)) and lim^ (x) = +oo. We have 

^"S(t)>l 
Integrating from to t, we obtain 

i/joS(t) >t,t> 

since ^ is decreasing 

S(t) < (t) ,t >0 

(2) 

(a) Let ip, the function defined by 

, , , /-P- Io r(o) ds 
ip(x)= / — — . 

Jx P (S) 

then ?/> is a strictly decreasing function on ]0,p _1 o T (0)1 and lim^ (x) = +oo. 
We have 

° rW = f^— 

from (|1.15j) , we infer that 

— ?/> o o r ct) > c 

(it 

Integrating from to t, we obtain 

^op- 1 oT(t) > ct 

this gives 

r (t) < p o V^ 1 (ct) , for every t > (1.17) 

Setting 

y (t) = KTp- 1 (Ct),t > 0. 

We have 

y' (f) + p (y (t)) = - CK p o t/;- 1 (ct) + p (k (^ _1 (ct))) 
Using (fTT3D and (fTTTTD 

?/' (t) + p (y (t)) > (mp (k) — ck) p o (ct) 
> (mp (k) — ck) T (t) 

(|1.16p gives 

y'(t)+p(y(t)) >T(t) 

y(o) >5(o) 

the result follows from the following lemma 

Lemma 1. Letpi (i = 1,2) a positive strictly increasing function onM+. Suppose 
that S and y are absolutely continuous functions and satisfy 

A Q 

"' +Pi(S) <T(t) on [0,+oo[. (1.18) 



dt 
and 

dy 



( , t +P2{y)>T 1 (t) on [0,+oo[. (1.19) 
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where r;Ti € ^([Ojoo)); I\ > T > 0; Pi > P2 > 0. in addition; if 

2/(0)>5(0) 

T/ien 

y (t) > S (*) ; /or t > 

First we finish the proof of the proposition, then we give the proof of the lemma, 
(b) Setting 

y (f) = Kp- 1 O r (t) , t > 0. 

We have 

y' (t) +p(v (t)) = k (p- 1 o r)' (t) + P (^p- 1 o r (t)) 

Using (|1.14p and the fact that 

^p" 1 (T (t)) + cF (t) > 
for some c > 0, we obtain 

y' (*) + p (y (*)) > {mp («) - ck) r (i) 

(|1.16p gives 

2/'(t)+p(y(t))>r(t) 
y(o)>s(o) 

the result follows from lemma [1] 

□ 



The proof of lemma [T] is borrowed from |27| 
Proof of lemmaUl Suppose that there exists to i n [0, +oo[ , such that 

S (to) = V (to) and S (t) > y (t) on [t , t + e] 
for some e > 0. Integrate ( 11.18P and (I1.19P from to to to + e > we obtain 



and 



therefore 



which gives 



S(t + e)-S(t )+ f° 6 pi(S(t))dt< f° C T(t)dt 

Jto Jto 

y(t + e)-y (t ) + I ° p 2 (y (t)) dt > f° Ti (t) dt 

Jto Jto 

S(t + e)+ I ° + p! (S (t)) dt<y(t + e)+ f ° + p 2 (y (t)) dt 

Jto Jto 



S(to + e)-y(to + e) < [° p 2 (y (t)) - Pl (S (t)) dt 

Jto 



i to 

< f° e pi (y(t))- Pl (S(t))dt<0 
Jto 



which contradict the fact that S (t) > y (t) on [to , to + e] • □ 
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Setting 

r (t) = 2 J m |/ (t, x)\ 2 dx + r(\\f (*, •) \\ lH m)) 

where ip* is the convex conjugate of the function ■0, defined by 
and 

^* ( s ) = sup [sy - 99 (y)] 



1 



1 



Superlinear damping: Assume 

s 2 e~^ < s < 1 
-s 2 e"^ -1 < s < 

We choose /iq 1 (s) = s 3 / 2 e~s, < s < rj « 1 and 

K >> max (E u (0) + ||r|| Ll(R+) , C T 

We have 

L 2 (M) 

The ODE (jl.5p governing the energy bound reduces to 

^ + C 75 3 / 2 e-^ <T(t) 
at 

with C > depends on E u (0) . If there are constants M > and 6* > 1, such that 

r(t) < M{\ + ty e 

then 

with c, Co, ci > 0. These constants may depend on E u (0) . 
Sublinear near the origin: Assume g (s) s ~ |s| 1+r ° , |s| < 1, ro € (0, 1) . We choose 
/i ( s ) = s 2 n>/(i+r ) for < s < 1 and 



» max (£„ (0) + ||r|| L i (R+) , C T ) ■ 

We have 

r (ii/ (*. •) ii^(m)) < c (ii/ (*, •) iistio + ( *' -)H' 2 w 

The ODE f| 1 . 5 1) governing the energy bound reduces to 

^ + C5 (l+ro)/2r < p (t) 

with C > depends on E u (0) . 

(1) If there are constants M > and > 1, such that 

r(t) < m(i + 

Then 
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(a) 9 € 



i l+rg 
' l-r 



E u (t) <c{l + t-T) ,t>T 



(b) > 1+m 



where c > 0. 

lo 
ro 

< c(t-T)" 1 ^ ,t > T 

with c > and depends on E u (0) . 
(2) If there are constants M > and 9 > 0, such that 

r (t) < Me _fl * 

Then 

E u (t) <c(t-Ty^ ,t> T 
with c > and depends on E u (0) 

2. The linear case: Proof of theorem [1] 
2.1. Preliminary results. 

Proposition 2. Zei u be a solution of (jl.ip mf/i initial data in the energy space. Then 

E u {t)< + ^ e^^Euis) + £ Jjf (a,x)\ 2 dxdo^j (2.1) 

/or every e > and /or every t > s > 0. 
Proof. Let £ > s > 0. From the energy identity 

E u (t) <E u (s)+ [ [ fd t udxda 

Js J M 

Using Young's inequality 

E u (t) < E u (s) + - / f |/ (a, x)\ 2 dxda + e f E u (a) da 

e Js JM Js 

for every e > 0. Now GronwalPs lemma, gives 

E u (t) < e< 1 -^ (e u (s) + IJ J \f (a, x)\ 2 dxda) . 

□ 

The result below is a generalisation of the comparison lemma of Lasiecka and Tataru [15j . 

Lemma 2. Let T > and 

• T € Lj oc (R+) and, non negative. Setting 5 (t) = f* +T T (s) ds. 

• W (t) be a non negative, continuous function for t £ R+. Moreover we assume that 
there exists a positive, monotone, increasing function a such that 



W(t) <a(t-s) 



W(s) + J r O) da 



for every t > s > 0. 



ENERGY DECAY 11 

• Suppose that £ and I — I: M + — > R are increasing functions with £(0) = and 

W ((m + 1) T) + £ {W (mT) + 5 (mT)} < W (mT) + 5 (mT) (2.2) 
for m = 0, 1, 2, .. where £ (s) does not depend on m. Then 

W(t) <a(T)(s(t-T) + 2 T{s)ds^j, Mt > T 

where S (t) is a positive solution of the following nonlinear differential equation 

^ + ±£(S) = T(t) ; S(0) = W(P). (2.3) 

Proof. To prove this result we use induction. Assume that W (mT) < S (mT) and prove that 
W ((m + 1) T) < S ((m + 1) T) where S (t) is the solution of (H3D • 
Integrating the equation (I2.3j) from mT to (m + 1) T yields 

i f (m+l)T 

S ((m + 1) T) = S (mT) - - / £(S (t)) dt + 5 (mT) (2.4) 

T JmT 



On the other hand, we have 

d_ / /■'.. , \ 1 

dt 

therefore, for t± > ti 



(s-J T( S )ds^=-^£(S)<0. 
S(h)<S(t 2 )+ f tl T(s)ds 



the function £ is increasing 

£(S(t)) < t(s(mt) + J T(s)ds^j , for mT < t < (m + 1)T 

< £ (S (mt) + 5 (mT)) 
Using now, (|2.4p , we obtain 

S ((m + 1) T) > S (mT) + 5 (mT) -£(S (mt) + 5 (mT)) 
Since the function, I — £ is increasing 

S ((m + 1) T) > W (mT) + 5 (mT) -£(W (mt) + 5 (mT)) 

$ZM , gives 

S((m + l)T) > W((m + 1)T). 
Setting t = mT + t, with < r < T. Then we obtain 

W(t) < a(r) 



W(t-r)+ [ F(s)ds 

Jt-T 

< a(r)(s(t-T) + J r(s)ds\ 

< a (T) (^S (t - T) + 2 jf" T (s) ds^j , for every t > T. 



□ 
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Proposition 3. We assume that (oj,T) satisfies the assumption (G). Then there exists Ct > 
0, such that the following inequality 



E u (t) < C T 



t+T 



t JM 



a (x) \dtu\ 2 + |/ (s, x)\ 2 alxds 



(2.5) 



holds for every t > 0, for every solution u of (jj.jp with initial data in the energy space *H, 
for every f in L 2 (R + x M) . 

Proof. To prove this result we argue by contradiction. We assume that there exist a sequence 
(u n ) n solution of (I1.7P with initial data in the energy space, a non- negative sequence (t n ) n 
and f n in I? (R + x M) , such that 



E Un (t n ) > n Jlf T f M a (x) \d t u n \ 2 + \f n (t, x)\ 2 dxdt, 
Moreover, u n has the following regularity 

u n G C (R+,H% (M)) n C 1 (R+,L 2 (M)) 
Setting a n = (E Un (tn)) 1 ^ 2 > and v n (t, x) = ""(^ +t ' rz ) , Then v n satisfies 

dfv n - Av n + a (x) d t v n = ±f n (t n + t,x) M + x M 
v n = K + x <9M 

(v„ (0) , <9 s u n (0)) = (« n (i n ) , (i n )) 



(2.6) 
(2.7) 

(2.8) 



Moreover, 



and 



(0) = 1 



1 > n J T J M a (x) |<9^ n | 2 + -±f n (t n + t, x) 
From the inequality above, we infer that 



dxdt, 



In lM a ( x )\dtv n \ dxdt -> 0, 

2 



Jo /a/ oL ^ n + X ; 



ctodi -> 

n— >oo 



We have 



Therefore, 



v n G C ([0, T] , Hi (M)) n C 1 ([0, T] , L 2 (M)) 



(2.9) 



(2.10) 



(2.11) 



E Vn {t) = E Vn (0)- [ [ a{x)\d t v n \ 2 dxdt+ [ [ —f n (t n + t,x)d t v n dxdt (2.12) 

JO JM JO JM a n 

and using (|2.ip , we infer that 

< 2e T (^„(0)+ f [ 

\ J0 JM 



'fn (tn ~t~ t, X i 



dxdt 



E Vn (t) < 2c 1 ( 1 + - ) , for every t G [0, T] 
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This estimate allows one to show that the sequence (v n ,d t v n ) is bounded in L°° ((0, T) ,H) 
then it admits a subsequence still denoted by (v n ,dtv n ) that converges weakly-* to (v,dtv) 
in L°° ((0, T) , %) . Passing to the limit in the system satisfied by v n , we obtain 

dfv-Av = ]0,T[xM 

d t v = ]0,T[xu (2.13) 

(v (0) , d s v (0)) G H 

and the solution v is in the class 

C ([0, T] , H l Q (M)) n C 1 ([0, T] , L 2 (M)) 

We deduce as in J. Rauch and M. Taylor [22] or C. Bardos, G. Lebeau, J. Rauch [1] that 
the set of such solutions is finite dimensional and admits an eigenvector v for A. By unique 
continuation for second order elliptic operator, we get dtv = 0. Multiplying the equation 
by v and integrating, we obtain v = 0. Now we prove that v n — > 0, in the strong topology 
of H~i oc ((0, T) (M)) . For that we use the notion of microlocal defect measures. These 
measures were introduced by P Gerard [TO] and L. Tartar [25]. Let \i the microlocal defect 
measure associated to the sequence (v n ) . From (I2.10p we infer that the supprot of fx is 
contained in characteristic set of the wave operator and it propagates along the geodesic flow 
(G. Lebeau [18]). Therefore 

Now the assumption (G) combined with the propagation of \x along geodesic flow, gives 

v n ^0, Hi oc ((0, T) , H 1 (M)) 

This gives 

/ <p(t)E Vn (t)dt -> (2.14) 
Jo 

for every (p 6 Cq 30 ([0,T]) . On the other hand, E Vn (0) = 1, therefore, from f|2. 12[) and the 
fact that 

(-T 



j j \d t v n \ 2 dxdt < 2Te T (l + ^) 



we deduce that E v (t) — > 1, for every t € [0, Tl . Since E Vn (t) < 2, by Lebesgue's dominated 

n— >oo 

convergence theorem 

r-T pT 

if it) E Vn (t) dt ->■ (f (t) dt. 







for every ip € Cq° ([0, T]) . We obtain a contradiction by choosing (p such that J Q T ip (t) dt > 
0. □ 

2.2. Proof of Theorem [JJ Let u be a solution of (jl.ip with initial data in the energy space. 
From the energy identity we have 

rt+T p rt+T p 

I I a (x) \d t u\ 2 dxdt = E u (t) - E u (t + T) +/ / / (s , x) d t udxds 
Jt Jm Jt Jm 
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therefore, using Young's inequality 

ct+T 



/ / a (x) \d t u\ 2 dxdt < E u (t) — E u (t + T) 
Jt Jm 



+e 



t+T 



M 



2 1 

\dtu\ dxds H — 



t+T 



\f(s,x)\ dxds 



M 



for every e > 0. Now using the observability estimate (|2,5p and (|2.ip , we can show that 



t+T 



M 



\d t u\ 2 dxds < 2Te T C T 



t+T 



M 



a (x) \dtu\ 2 + \f (s, x)\ 2 dxds 



(2.15) 



Then setting e 



4Te T C T 



we infer that 



t+T 



M 



a (x) \d t u\ 2 dxdt < E u (t) - E u (t + T) + ( 1 + Te T C T 



t+T 



A I 



\f(s,x)\ dxds 



Hence with C T = ( 1 + Te 1 C T ) , we have 

ct+T 



a{x) \d t u\ z dxdt < 2 



M 



E u (t) -E u (t + T) + C T 



t+T 



|/ (s, x)\ 2 dxds 



M 



Now from the observability estimate (|2.15p 



E u (t)-E u (t + T) + (C T + 1 



E u (t) < 2C T 



with C T > 1. Setting C 1)T = 2 ( C T + 1J . We remark that 

ft+T 



t+T 



\f(s,x)\ dxds 



M 



2C T 



E u (t) -E u (t + T) + C x ,t 



M 



\f(s,x)\ dxds 



> C x 



t+T 



M 



\f(s,x)\ dxds 



Then for Ct = ^Ct, we have 

i-t+T 



K(t) + Ci, T [ [ \f(s,x)\ 2 dxds < C T [E u (t)-E u (t + T) 
Jt Jm 



t+T 



\f(s,x)\ dxds 



M 



Therefore 



E u (t + T) + — 



E u (t) + C X ,T f +T [ \f(s,x)\ 2 dxds 
Jt Jm 

[■t+T r 

< E u (t) + C hT / / \f(s,x)\ 2 dxds 
Jt Jm 

Setting t = mT, with m € N and using Lemma |2, 41 we conclude that 

E u (t) < 4e T (s(t-T) + £ T (s) ds^j , Vt > T 
where S (t) is a positive solution of the following nonlinear differential equation 

^ + -^fS = T(t) ; S(0) = E u (0) 



(2.16) 
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and 

r( s ) = c 1)T / \f(s,x)\ 2 dx 

Jm 

3. The nonlinear case: Proof of theorem [2] 

This part is devoted to the proof of theorem [2j First we give the following energy inequality. 

Proposition 4. Let u be a solution of ( | with initial data in the energy space. Then the 
following inequality 

E u {t)< (l + Vje^ (e u (s) + £ Jjfi^xrfdxdo-^ (3.1) 

holds for every e > and for every t > s > 0. 



For the proof of (13.ip , we have only to proceed as in the proof of (12. ip . Now we give the 
proof of theorem [2j 

Proof of Theorem^ Let u be the solution of (jl.7p with initial condition (no, u\) in the energy 
space %. Let t > and <fi = u (t + •) be the solution of 

d 2 s (j) -A<j) + a(x)g {d s <t>) = f(t + s,x) R+ x O 

= M.+ x dn (3.2) 

(0(O),a s 0(O)) = (n(t),a t n(t)) 

We argue as in [9]. Define z = <f> — v, where v is the solution of (jl.lj) with initial data 
(t) , 5 t n (t)) and / = /(* + •,•)• Then z satisfies the system 

' <9 t 2 z - Az + a (x) g (d t <p) - a(x)d t v = R + x!l 
< z = R + xffl 
k (z(0),d t z(0)) = 

Let T > 0, such that (to, T) satisfies the assumption (G) . It is clear that a (x) [g {dt4>) — dtv) € 
L 2 ((0,T) x Q) . This observation permits us to apply energy identity, whence 

Ez(T) = [ [ a(x){d t v-g{d t <j)))d t zdxdt 

a (x) ( \d t v\ 2 + g{dt4>)dt4>) dxdt + g (d t 4>) d t v + d t vd t <pdm a 

o Jm v 7 Jo Jm 

The monotonicity of g (g (s) s > 0) and the estimate above, gives the following estimate: 

/ / a{x)(\d t v\ 2 + g{d t <j>)d t (p) dxdt < \ \ g (d t (j)) d t v + d t vd t (j)dm a 
Jo Jm v y Jo Jm 

Now the observability estimate (|2.5p . gives 

E u (t) = E v (0) < C T ( C I 9 {d t <j>) dtv + d t vd t 4>dm a + T I \f(t + s,x)\ 2 dxds 
\Jo Jm Jo Jm 

for some Ct > 1. From the estimate above we infer that 

E u (t)<C T ( [ + [ g (d t u) 8 t v + d t vd t udm a + f + [ \f (s,x)\ 2 dxds) 
\Jt Jm Jt Jm J 

where v (s) = v (s — t) , s > t > 0. 



o Jm 
t 
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Lemma 3. Setting 

M a>t = [s,t]xQ, t>s>0 and M 0jt = M t 
Let t>0. For i = 0,1 let 

M° = {(s, x) E[t,t + T]x ft; \d s u (s, x)\ < Vo } , M\ = M t>t+T \ M° 

and define 

6 (M tit+t ) = [ \d s u (s) d s v (s - t)\dm a , (M*) = [ \g (d s u (s)) d s v (s - t)\ dm ai 

JM t J Ml 

where u and v denotes respectively the solution of 7[ ) and (jl.ip with initial data (uq,ui) 
and ((vo,vi) = u(t) ,d t u(t)). 

(1) The following inequality holds for every e > 0, 

* (M f °) + 6 (M°) < eE u (t) + -m Q (M T ) ( ^ f 9 (d s u) 8 s u dm a ) (3.3) 

e \m a {M T ) J Mt t+T J 

rt+T r 

+e / / \f(s,x)\ 2 dxds (3.4) 
Jt Jm 

with C > 0. 

(2) Estimate on the damping near infinity. The following inequality 

* (Ml) + 6 (Ml) < eE u (t) + Ce- 1 ( / g (d s u) 8 s u dm a ) (3.5) 

yM t:t+T J 



pt+T 

+e / / \f(s,x)\ 2 dxds (3.6) 
Jt Jm 

holds for every e > with C > 0. 

For the proof of the lemma above we have only to proceed as in [9j Lemma 3.1 cases 1 and 
2] and to use the energy inequality (|2.ip . 
Now using (|3.3p and (|3.5p , we deduce that 

E u (t) < C T [ eE u (t) + C T , e h I [ g (d s u) d s udm a + f \f (s, x)\ 2 dxds ) ) 

\ \JMt t t+T JM t<t +T J J 

for every e > 0, where the function h = I + m a (Mt) ho o — and Ct > 1- Setting e small 
enough, e.g. e = ^ 

E u (t)<C T h( [ + [ g(d s u)d s udm a + [ + [ \f (s,x)\ 2 dxds) (3.7) 
\Jt Jm Jt Jm J 

for some Ct > 1. This gives 

E u {t)+ / \f(s,x)\ 2 dxds<2C T h( / g (d s u) d s udm a + / [/ (s, x)| 2 dxds 

(3.8) 

On the other hand, the energy identity gives 

rt+T /• rt+T r 

\ \ a (x) g (d t u) d t udxda < E u (t) - E u (t + T) + / |/ (a, x) d t u\ dxda (3.9) 

Jt Jm Jt Jm 
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Let ifi, defined by 

{ +00 s G M*_ 

It is clear that ip convex and proper function. Hence, we can apply Young's inequality |23j 

t+T r rt+T 

\f (a,x)d t u\dxda < / ||/ (a, .) \\ L2 \\d t u (a, .) \\ L2 da 

M Jt 

< f t + r(\\f(.v,.)\\ L i) + 4>(\\dtu(a,.)\\ L2 )da 

where ip* is the convex conjugate of the function tjj, defined by 

^* ( s ) = S up [sy - ip (y)] 



Using the energy inequality (|3. 1 1) and the observability estimate (|3.8[) , we infer that 

't+T ^ / /-t+T r rt+T p \ 

^(||9t«(cr, .)\\ L 2)da < - I / / g (d s u) d s udm a + / / |/ (s, x)| 2 dxds ) 

2 \Jt jm Jt Jm J 



then (|3.9p , gives 

rt+T r / rt+T r 

/ / a (x) g {d t u) d t udxda < 2 ( E u (t) — E u (t + T) + / / \f (s, x)\ 2 dxp.W) 
Jt Jm V Jt Jm 

+ r{\\f{^.)\\ L ,)da^ (3.11) 

The inequality above combined with the observability estimate (|3.7|) and the fact h = 
I + m a (M T ) h o ma ( MT ) is icreasing, gives 

E u (t) < C T h ( 4 ( E u (t) - E u (t + T) + f +T [ \f(s,x)\ 2 dxda + ^ (||/ (a, .)\\ L2 ) da 

V V Jt Jm Jt 

Setting 

r» = 2 [ \f( S ,x)\ 2 dx + r(\\f(s,.)\\ L 2) 
Jm 

Therefore 

E u (t) + £ + r (s) ds < Kh U (e u (t)-E u (t + T) + J + T (s) ctecfej 



with K>C T . Setting (t) = f/ +T T (s) (is. Thus 

E u (t + T) + hi' 1 (± (E u (t) + 9 (t))) < E u (t) + 9 (t) 
for every t > 0. Take i = mi, m £ N 

£ u ((m + 1) T) + i/r 1 f -I (£„ (mT) + 9 (mT))j < (mT) + (mT) 
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Setting W (t) = E u (t) ,£(s) = ^ o ^ and T (t) = 2 J M \f (s, x)\ 2 dx + ^dlf (*, •) M • It 
is clear that the functions £ and I — £ are increasing on the positive axis and I (0) = 0. The 
function T € Lj and non negative on M + . According to lemma [2] 

E u (t) < 4e T (s(t-T) + J r (s) ds^ , Vt > T 
where 5 (t) is the solution of the following nonlinear differential equation 

^ + ±£(S)=r(t) ; S(0) = W(0). (3.12) 

□ 
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